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Abstract 

We show that a high frequency standing wave in SU(2) gauge theory is 
unstable against decay into long wavelength modes. This provides a non- 
perturbative mechanism for energy transfer from initial high momentum 
modes to final states with low momentum excitations. The Abelian case does 
not manifest such instability. Our analysis is supported by lattice simulations. 
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Recently the intriguing problem of the possibly large rate for electroweak baryon-number 
changing processes in high-energy collisions, originally based on the instanton approach [|TJ , 
was considered from a completely different point of view, based on the classical analogue to 
high-energy particle collisions 0|| . The hope was that one could observe the energy transfer 
from fast (high frequency) classical wave modes, presumably corresponding to initial states 
containing few high-energy particles, to slow (low frequency) classical wave modes, repre- 
senting low-energy multi-particle final states. The studies involved numerical simulations 
of the (l+l)-dimensional abelian Higgs model [[J and the 4 -theory in (3+1) space-time 
dimensions [Q]. In both cases no indication for the existence of such a mechanism was found. 
Instead, the classical high-energy scattering behavior was found to be fully consistent with 
expectations from lowest order perturbation theory. The nonlinearity of the investigated the- 
ories apparently does not furnish a mechanism for the formation of final states containing 
many low-energy particles from the initial states with few high-energy particles. 

In these approaches the main object of study was the nature of the transfer of energy 
from high-energy to low-energy modes, which is thought to be in some sense analogous to 
the high-energy scattering situation. We believe that the observed absence of an efficient 
energy transfer mechanism is intimately connected with the integrable nature of the consid- 
ered classical systems. On the other hand, it is well known that nonabelian gauge theories 
are nonintegrable in the classical limit [[§-§]. So in these theories there might exist nonper- 
turbative classical mechanisms for the coupling between fast and slow modes. In this paper 
we want to investigate the existence of such a mechanism in the framework of a simple 
instructive example involving a nonabelian gauge theory. 

Let us model a high-energy particle collision as the interaction between two counter- 
propagating Yang-Mills plane waves. For simplicity let us work in SU(2) gauge theory. The 
simplest ansatz satisfying the free Yang-Mills equation is an abelian standing wave 

A1(x,t) = S&S&A cos k x cosc^ot (1) 

with k = ujq assumed to be large. (We work in the temporal gauge Aq = 0.) We now consider 
small perturbations a^(x,t) around this solution and study their stability properties. We 
will show that there are unstable low-energy modes for arbitrarily small amplitude A of the 
driving wave. The amplitude of these modes grows exponentially with time with a growth 
rate proportional to A. We will later confirm the presence of the instability by a numerical 
calculation and demonstrate that it leads to rapid total thermalization of the field energy, 
once the perturbation has grown to the same strength as the driving wave. 

It is easy to see that fluctuations in the same color direction as the background field 
(c = 3) are stable, hence we will consider only color fluctuations transverse in color space 
a, L = a} + iaf. The linearized field equations and Gauss' law for the field perturbation are 

d 2 t a = (D ■ D)a - 2i(B x a) - D(D ■ a) (2) 
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and 



d t (D ■ a) = 0, (3) 

where boldface type indicates three dimensional vectors in position space. B is the mag- 
netic background field generated by the standing wave, and D = V + iA is the gauge 
covariant derivative. Gauss' law (3) ensures that the physical states are invariant under 
time-independent gauge transformations. To fix the residual time-independent gauge, we 
impose the background field Coulomb gauge constraint: 

D ■ a = 0. (4) 

Now the field equation simplifies, allowing for a decomposition of the orientations of the 
polarization vector of the perturbation. Modes polarized in the direction of the magnetic 
background field B (the y-direction) are stable; hence we only need to consider modes of 
the form 

a(x, t) — (z + iax)a a (x, t), (a = ±1) (5) 

where x and z denote unit vectors in the direction of the wave vector and the polarization 
of the background field, respectively. These modes obey the equation 

d 2 a a - V 2 a ff = (2iA ■ V - A 2 - 2aB)a a . (6) 

We now note that the background field A\ oscillates very rapidly in space and time, 
whereas we are seeking slowly varying perturbative modes a\. The dynamics of these modes 
is governed by the time- and space-averaged interaction with the fast background field. Since 
we consider the amplitude of the background wave to be small, i.e. A <C ko, we will find a 
solution perturbatively by calculating the selfenergy of the fluctuating field to second order 
in the background field. After spacetime averaging, we obtain the following expression for 
the selfenergy for a fluctuation with wave vector |fe| = k ko and frequency u ujq\ 
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, 2 u 2 + kl 



k z ) 



(7) 



Note that the expression for E is independent of the sign of the circular polarization a of 
the perturbation ([|). The dispersion relation for the fluctuation modes is obtained from the 
poles of the background field propagator 

G(k, uj)- 1 = -u 2 + k 2 - E(fe, u). (8) 

We shall check that the perturbative solutions also satisfy the background Coulomb 
gauge constraint (0). Since the equation of motion M) is solved to second order in A, it is 
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consistent that we only require the solutions to satisfy Gauss' law to second order. For a 
particular momentum k, we look for a solution of the form 



a(x,t) = (z + i<jx)a a {x, t). (9) 



<T=±1 



Inserting it into (f|), disregarding the fast oscillating terms, we get a relation between 
and a_, which is represented by their ratio as a function of momentum k, 



where 




(10) 



Hence for any k there is a solution satisfying Gauss' law. The most unstable modes are 
those with k y = k z = for which, keeping only terms up to second order in A, we obtain 
from (H) the following dispersion relation: 

u 2 = kl±-j=Ak Xi (12) 

which is complex for any non-zero k x and A. We conclude that the infrared instability 
exists for arbitrarily small values of the amplitude A of the background standing wave. 
When k x 3> A, the imaginary part of the frequency, i.e. the exponential growth rate of a 
perturbation, 

M"(*)l = (13) 

is independent of the wave vector k x . 

Unfortunately, the convergence of the above expansion in increasing powers of A 2 is not 
evident. An explicit calculation shows that the fourth order contribution to selfenergy has 
the same form as the second order one (7) near the pole, but with a smaller numerical 
factor. In order to verify the existence of this instability beyond perturbation theory, and 
to find out what happens when the fluctuation begins to absorb a significant fraction of 
the driving background wave, we have studied the evolution of a slightly perturbed low- 
amplitude standing plane wave in SU(2) gauge theory on a three-dimensional lattice in the 
classical limit. The numerical aspects of such simulations have been described in detail 
elsewhere ||[7|. Here we have initialized the gauge field as an abelian standing wave with 
an amplitude corresponding to less than 5% of the maximal magnetic energy density on the 
16 3 lattice. When we apply a small abelian perturbation, restricted to the same direction in 
color space as the standing wave, the field oscillations remain stable, as shown by the dashed 
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line in Fig. 1. When we add a general nonabelian perturbation, pointing in a random color 
direction, the field oscillations develop a visible instability around time t = 20, as shown by 
the solid curve in Fig. 1. 

It is instructive to compare the Fourier spectrum of the magnetic field energy density 

E m (k) = ~J B(x) 2 e* kx d 3 x, (14) 

at the end of our simulation (t = 50) with the energy spectrum of the initial configuration, 
as shown in Fig. 2, in both cases. The initial spectrum is almost completely concentrated 
in the mode with \k\ = it /a, where a is the lattice spacing. The case with a non-abelian 
perturbation is shown in the upper part of Fig. 2, where the final spectrum is distributed 
over all modes. On the other hand, in the case with an abelian perturbation, the final 
distribution remains concentrated in the original mode, as shown in the lower part of Fig. 
2. In order to check whether the available field energy has been thermalized in the case with 
non-abelian perturbation, we have calculated the probability distribution of the magnetic 
plaquette energy P(E), divided by the single plaquette phase space u(E). The distribution 
P(E)/ar(E) falls exponentially with E, indicating thermalization Jf|. 

The quantitative feature of the instability can be characterized by the associated largest 
Lyapunov exponent Ao- Numerically we find, in dimensionless form, 

A a = kaf(a), a = — — , (15) 

ka 

where g is the coupling constant. Obviously, Ao survives in the continuum limit a — > 0. The 
function / is obtained numerically and shown in Fig. 3. For small a it is a linear function 
/(en) ~ 0.5a. Inserting into the above relation, we have 

Ao « (16) 



which resembles our result (13|) from second order perturbation but involves a different 
coefficient. We note that the scaling behavior fll~5|) is different from the scaling behavior of 
a Lyapunov exponent of a random trajectory [[7], which scales with the energy rather than 
the amplitude of the background field. The unstable mode described above is more directly 
related to the instability of a uniform chromomagnetic background field f8[, and can be 
traced back to the magnetic moment interaction of a spin-1 particle, expressed by the form 
2a B in (|). 

Let us make one final remark connected with the high-energy limit of eq. @ for the 
perturbation a(x,t). If we consider the behavior of this equation under rescaling of the 
longitudinal coordinates x and t (or that of the light-cone coordinates x ± t) by 

t — > At, x — * Xx, y y, z — > z, (17) 
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one will see that after taking the high-energy limit || A — > the terms with magnetic field 
B in (H), which caused the above described instability, will disappear. One is left with the 
trivial wave equation 

- flg)a(iM) = 0. (18) 

This example shows that the transverse directions must be handled with care because some- 
times they can be the origins of singularities associated with the gauge sector of nonabelian 
gauge theories. In our example the color field does not behave like an abelian electromagnetic 
field, but shows the complexity of the transverse dynamics that cannot be ignored. 

In conclusion, we have shown that a standing wave of high frequency in SU(2) gauge 
theory is unstable against decay into long wavelength modes. This implies that interactions 
between high energy particles can result in a final state with many low energy particles, 
which might be relevant in baryon-number changing process. 
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Figure Captions 



Fig.l: The total electric energy as a function of time. The initial state is a standing wave 
with a small perturbation. The solid line shows the case when the perturbation is non- 
abelian, which results in the destruction of the plane wave beginning at t ~ 20. The 
dotted line is for an abelian perturbation, in which case the standing wave is stable. 
The computation is performed on a 16 3 lattice and energy per plaquette is about 0.16 
in lattice units. 

Fig.2: Time evolution of energy spectra. The upper part is for the case of a non-abelian 
perturbation, against which the original standing wave is unstable, while the lower 
part is for an abelian perturbation, where we observe no instability. 

Fig.3: The scaling function f(ct) defined in eq.(|l5|). The dotted line is a linear fit with slope 
0.5. The results are obtained by varying the amplitude A with a fixed wave number 
ka = 0.49 on a 256 x 4 2 lattice. 
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